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model  provided  attenuation  results  in  close  agreement  with  those  computed 
using  full  Mie  theory,  whila  avoiding  the  computational  complexities  associated 
with  a detailed  treatment  of  scattering  effects. 

With  the  Gross-Reber  line  shape  ~nd  the  precipitation  model  incorporated 
in  the  simulation  programs,  sets  of  channels  in  two  60  GHz  oxygen  complex 
and  near  the  118  GHz  oxygen  line  were  selected,  simulated,  and  inverted  for 
atmospheric  conditions  representing  mid  latitude  spring  and  tropical  summer. 

The  results  of  the  inversions  from  each  set  of  channels  were  quite  good 
through  the  118  GHZ  were  slightly  less  accurate  at  low  altitudes  due  to  a 
greater  sensitivity  to  water  vapor.  Separate  inversion  analyses  for  specific 
cloud  and  h/drometeor  conditions  showed  that  the  60  GHz  channels  were  more 
accurate  for  temperature  profiling  in  the  presence  of  clouds  and  that  temper- 
ature inversion  results  in  the  presence  of  rain  were  seriously  degraded  for 
both  sets  of  channels. 

Inversion  for  temperature  data  using  these  channels  in  the  uplooking 
mode  showed  excellent  results  up  to  500  mb,  without  any  channel  optimization. 
Separate  investigations  of  the  use  of  these  oxygen  channels  with  channels 
near  the  22  GHz  and  183  GHz  water  vapor  lines  showed  great  promise  in  tht. 
determination  of  cloud  top  height,  cloud  phase,  layered  liquid  water  and 
prop-size  distribution. 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGEOWi en  D«l»  Entered) 


ACKNOWLEDGMENTS 


The  authors  wish  to  acknowledge  the  valuable  assistance  provided  by 
Mr.  James  H.  Willand  and  Ms.  Nancy  Tripp  in  the  computer  modeling  and  in 
the  simulations  experiments.  They  also  wish  to  thank  Ms.  Cha-Nan  Tang 
for  her  services  in  data  reduction  and  analysis. 


TABLE  0!-'  CONTENTS 


I’age 


A C K N 0 W l . I . I > CM  1 . N ’ l S 

1 . INTRODUCTION 

1 . 1 Purpose 

1 . 2 Scope 

2.  PRECIPITATION  MODEL 

2.1  Background 

2.2  Theoretical 

2.3  Rayleigh  Limit  - Stael^n  Approximation 

2.4  The  ERE  Precipitation  Extinction  Approximation 

2.5  Drop  Size  Distribution  lunction 

2.6  Discussion  of  the  Approximation 

3.  TEMPERATURE  PROFILING  ANALYSIS 

3 . 1 Scope 

3.2  Approach 

3.2.1  Line  Shape  Analysis 

3.2.2  Temperature  Profiling  Analysis 

3.3  The  Attenuation  Due  to  Oxygen 

3.3.1  The  Oxygen  Line  Shapes 

3.3.2  A Comparison  Between  the  Van  Vleck-Weisskopf 
and  the  Gross  Line  Shapes 

3.4  Inversion  Analysis  Results 

3.4.1  Selection  of  Spectral  Channels 

3.4.2  Temperature  Inversion  (vaults:  The  General 

Case 

3.4  3 Effects  oi  _ouds  and  Rain  on  Temperature 
P rofiling 

3.5  Effects  of  Droplet  Size 

4.  PROBING  EOR  CLOUD  PARAMETERS 

4.1  I'revious  Results 

4.2  Liquid  Water  Profile 

4.3  Cloud  Type  Discrimination  and  Cloud  Top  Altitude 
Determination 

5.  SUMMARY  AND  CONCLUSIONS 

6.  REFERENCES 


1 1 
1.3 
ID 
19 
19 
19 
19 
2 n 
> -> 


31 


55 

38 
6 2 
64 


v 


Preceding  page  blank 


LIST  or  F I CURIES 


2-1  Range  of  Prop  Sizes  and  Mie  Hfficiency  factors  at 

X = 1 cm  1 •’ 

2-2  Reduction  of  the  LRT  Approximation  to  the  Staelin 

Formula  as  the  Rainfall  Rate  Approaches  Zero  IS 

2-3  Comparison  of  the  [Extinction  Coefficient  Using  the 

Staelin  Formula,  F:RT  Approximation  and  Mie  Theory  16 

2-4  Comparison  Between  Mic  Theory  and  the  F.RT 

Approximation  l(l 

2-5  Frequency  Dependence  of  the  [Extinction  Coefficient 

for  Various  Models  1 


3-1  flow  Diagram  for  Application  of  the  Data  Inversion 

Methodology  to  Parameter  [Elimination  for  a Radiometric 
F.xperiment 

3-2  Atmospheric  Oxygen  Attenuation,  Van  Vleck-Wei sskopf 
Line  Shape 

3-3  Atmospheric  Oxygen  Attenuation,  Gross  Line  Shape 

3-4  Comparison  of  the  Van  Vleck-Wei sskopf  and  Gross 

Temperature  Weighting  Functions  for  the  53.29  GHz 
Channel 


3-5  Comparison  of  the  Van  Vleck-Wei sskopf  and  Gross 
Temperature  Weighting  Functions  for  the  54.3  Gllz 
Channel 


3-6  Temperature  Weighting  Functions  for  the  Channels  in 
the  60  GHz  Complex  Computed  for  Midlatitude  and 
Tropical  Atmospheres 

3-7  Temperature  Weighting  Functions  for  the  Channels 

Near  the  118  GHz  Line  Computed  for  Midlatitude  and 
Tropical  Atmospheres 

3-8  Temperature  Weighting  Functions  for  the  Channels 
Near  the  118  GHz  Line  Computed  for  a Dry  Tropical 
Atmosphere 

3-9  Inversion  Results  for  the  Midlatitude  Case  With 
General  Cloud  Statistics  Using  Channels  in  the 
60  GHz  Complex 


20 

25 

25 


26 


26 


28 


29 


30 


o.i 


vi 


LIST  or  FI CURLS  (Continued) 


•'igure  Page 

3-10  Inversion  Results  for  the  Midlatitude  Case  with 
General  Cloud  Statistics  Using  Channels  Near  the 
118  GHz  Line  34 

3-11  Inversion  Results  for  the  Tropical  Case  with 

General  Cloud  Statistics  Using  Channels  in  the 

60  GHz  Complex  35 

3-12  Inversion  Results  for  the  Tropical  Case  with 

General  Cloud  Statistics  Using  Channels  Near  the 

118  GHz  Line  3o 

3-13  Inversion  Results  for  the  llplooking  Midlatitude 

Case  Using  Channels  in  .he  60  GHz  Complex  38 

3-14  Inversion  Results  for  .he  llplooking  Midlatitude 

Case  Using  Channels  N' ar  the  118  GHz  Line  39 

3-15  Temperature  Weighting  Functions  of  the  Low 

Altitude  Channels  in  the  Presence  of  a 4 Km  Cloud  4f) 

3-16  Temperature  Weighting  Functions  of  the  Low 

Altitude  Channels  in  the  Presence  of  an  8 Km  Cloud  4 1 

3-17  Temperature  Weight '.ng  Functions  of  the  Low 
Altitude  Channels  in  the  Presence  of  a 6 Km 
Glaciated  Cloud  4 3 

3-18  Temperature  Weighting  Functions  of  the  Low 

Altitude  Channel1,  in  the  Presence  of  a 12  Km 

Glaciated  Cloud  44 

3-19  Temperature  Weighting  Functions  of  the  Low 

Altitude  Channels  iri  the  Presen  •’  f mii  of 

Rain  4 5 

3-20  Temperature  Weight  1 ig  Functions  of  the  Low 

Altitude  Channels  in  the  Presence  of  4 Km  of 

Ice  and  8 Ijn  of  Rain  46 

3-21  Temperature  Inversion  Results  for  Clear  Sky, 

2 Km  Cloud,  and  2 Km  Rain  Cases  Using  Channels 

in  the  60  GHz  Complex  4S 

3-22  Temperature  Inversion  Results  for  2-4  Km  Cloud 
and  2-4  Km  Rain  Cases  Using  Channels  in  the 
60  GHz  Complex  4^ 


v i i 


LIST  OF  FIGURES  (Continued) 


Figure  Page 

3-23  Temperature  Inversion  Results  for  Clear  Sky, 

2 Ism  Cloud  and  2 Km  Rain  Cases  Using  Channels 

Near  the  118  GHz  Line  49 

3-24  Temperature  Inversion  Results  for  2-4  Km  Cloud 
and  2- 4 Km  Rain  Cases  Using  Channels  Near  the 
118  GHz  Line  49 

3- 25  The  Dependence  of  Total  Brightness  Temperature 

on  Cloud  Drop  Size  and  Microwave  Frequency 

4- 1  Microwave  and  Infrared  Measurements  of 

Glaciated  and  Ungiaciated  Clouds  °4 


vii  i 


LIST  OF  TABLES 


able  Page 

3- 1  Frequencies  Selected  for  Analysis  of  the  General 

Case  31 

4- 1  Inversion  Results  for  Integrated  Water  Vapor  and 

Layer  Integrated  Liquid  Water  57 

4-2  Summary  of  the  Temperature  Characteristics  of  Cloud 

Models  Used  for  Analysis  59 


1.  INTRODUCTION 


1 . 1 Purpose 

On  11  December  1972,  the  Nimbus  5 satellite  was  successfully  launched 
into  a near  circular  orbit  carrying  the  five-channel  Nimbus  I:  Microwave 
Spectrometer  (NKMS) . Since  launch,  this  radiometer  system  has  transmitted 
i large  volume  of  data  from  which  atmospheric  temperature  profiles  and 
integrated  water  vapor  values  have  been  successfully  inferred.  This  expe'i- 
ment , the  NEMS,  has  demonstrated  v.um.lusively  that  microwave  radiometry 
is  a viable  tool  to  obtain  global  m *asurenu:nts  of  atinospheri c parameter1' 
from  orbiting  spacecraft  on  a routine  basis. 

In  recognition  of  the  potentiils  of  such  a measurement  system  to  sat i f 
its  needs  for  global  atmospheric  lata,  the  li.  S.  Air  Force  initiated  a 
program  (under  Contract  No.  F047  1-72-C-0227)  to  define  the  configuration 
for  a system  which  can  best  provide  the  Air  Force  with  operational  data  on 
a global  basis,  within  the  constraints  of  weight  and  power  limitations 
imposed  by  the  spacecraft  used  in  the  Air  Force's  Defense  Meteorological 
Satellite  Program  (DMSP) . Environmental  Research  6 Technology,  Inc.  (ERT) 
contributed  to  this  program  by  performing  analytical  studies  of  the  expected 
performance  of  a system  configuration  under  a number  of  geophysical  condi- 
tions (Gaut,  Reifenstein,  Chang,  and  Blinn,  1973).  The  results  of  this 
study  demonstrate  the  usefulness  of  a microwave  system,  even  within  tl  ■ 
severe  limitations  placed  cn  the  system  design  by  the  spacecraf!  l ns. 
appears  that  the  Air  Force  will  proceed  with  the  procuremen’  actual 
microwave  radiometer  systems  for  launch  in  tb  ,r'"_  c.me  frame. 

The  system  under  consideration  1 .’  ..r  Force  is  based  on  the  state 

of  the  art  in  coaipor  :;;t  ..icuu.i  mg  techniques  and  on  the  current  analytical 

understanding  of  the  interaction  between  microwave  radiation  and  atmospheric 
parameters.  \ number  of  developments  have  taken  place  which  suggest  that  the 
concepts  in  microwave  remote  sensing  of  atmospheric  parameters  need  to  be 
further  analyzed  and  refined  to  advance  the  state  of  the  art  in  anticipation 
of  future  generations  of  sensors.  The  purpose  of  this  study  is  to  develon 
models,  analysis  techniques,  and  to  perform  analyses  of  a number  of  special 
problems  in  microwave  remote  probing  of  atmospheric  parameters.  The  results 
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arc  useful  in  interpreting  data  expected  from  the  first  microwave  sensor 
from  the  DMSP  and  in  designing  future  systems. 

1 . 2 Scope 

The  study  addresses  itself  to  two  main  topics  of  research  and  analysis: 
the  effects  of  clouds  and  hydromei  eors  on  temperature  profiling,  and  the 
possibility  of  obtaining  cloud  information  from  microwave  remote  sensing. 
However,  in  the  course  of  achieving  these  objectives,  other  developments  in 
CRT's  analytical  procedures  were  made.  The  most  significant  of  these 
i ncl  tide : 

• Implementation  of  rhe  Gross-Reber  line-shape  model  for  the 
oxygen  absorption  coefficients 

• hxtension  of  analysis  of  temperature  inversion  accuracies  to 
the  118  GHz  line 

• Development  of  an  approximate  model  for  the  analysis  of  the 
effects  of  precipitation 

e Analysis  of  the  application  of  the  183  GHz  line  shape  for 
atmospheric  probing 

The  results  of  these  efforts  are  organized  into  three  major  sections 
in  this  report.  In  Section  2,  the  precipitation  model  is  discussed.  The 
discussion  includes  the  assumptions  and  limitations  of  the  model.  Compu- 
tations made  us;'ig  the  model  are  compared  with  results  using  the  full  Mie 
approach  and  with  results  using  the  Staclin  approximation  for  the  analysis 
of  cloud  effects. 

Detailed  studies  of  the  problems  associated  with  temperature  inversions 
are  presented  in  Section  3.  The  differences  in  the  weighting  functions  and, 
therefore,  in  temperature  inversion,  are  first  discussed  for  selected 
frequencies.  Inversion  results  using  channels  in  the  60  GHz  complex  arc 
tiien  compared  with  results  obtained  for  the  118  GHz  line.  The  effects  of 
clouds  and  precipitation  on  temperature  sounding  in  each  of  the  two  oxygen 
absorption  regions  are  shown  by  actual  inversion  results  using  ensembles 
of  data,  with  and  without  hydrometeors,  finally,  the  effects  of  differing 
Crop  size  distributions  are  separately  analyzed  and  presented  using  the  full 
Mie  approach. 


In  Section  *1  the  resulto  of  analysis  of  remote  sensing  of  clouds  bv 
microwav<_  radiometry  are  presented.  Various  combinations  of  channels  were 
analyzed  in  an  attempt  to  define  the  parameters  which  could  be  inferred. 

The  parameters  included  liquid  cloud  top  height  and  ice  cloud  top  height, 
cloud  liquid  water  pr,file,  and  drop  size  parameters.  In  these  analyses, 
use  was  made  of  the  183  GHz  water  vapor  line. 

Many  of  the  problems  and  techniques  examined  in  this  study  have  never 
previously  been  explored  in  detail.  Thus,  the  results  provided  in  this  study 
are  preliminary.  Conclusions  reached  on  the  basis  of  these  results  are  sum- 
marized in  Section  5. 
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2.  PRECIPITATION  MODEL 


2.1  Background 

Precipitation  droplets  are  primary  contributors  to  nonresonant  opacity 
in  the  microwave  region  of  the  spectrum.  As  a consequence,  the  presence  of 
rain  may  severely  limit  the  accuracy  of  atmospheric  temperature  profile 
retrieval  from  microwave  measurements.  Evaluation  of  the  effects  of  precipi- 
tation on  temperature  measurements  by  "exact"  means  is  at  best  difficult  to 
perforin  due  to  the  computational  complexities  involved.  Therefore,  the 
puipose  of  the  work  described  in  this  section  is  the  development  oE  an 
approximation  to  the  interaction  which  would  permit  useful  analysis  of 
precipitation  effects. 

The  model  is  basically  an  analytical  extension  of  an  existing  approxi- 
mation model  for  the  analysis  of  clouds.  The  development  recognizes  that 
the  essential  difference  between  precipitating  and  nonprecipitating  clouds 
is  in  the  nature  of  their  respective  drop  size  distributions.  The  radiative 
properties  of  nonprecipitating  clouds  may  be  derived  by  assuming  individual 
cloud  droplets  to  be  Rayleigh  scatterers  and  integrating  over  the  range  of 
cloud  droplet  sizes.  This  leads  to  an  analytical  form  for  the  cloud  extinc- 
tion coefficient  whose  properties  in  the  centimeter  wavelength  region  have 
been  evaluated  (Staelin,  1966;  Gaut  and  Reifenstein,  1970).  Unfortunatel v, 
precipitation  droplets  are  one  to  two  orders  of  magnitude  larger  than  cloud 
droplets,  and  the  Rayleigh  approximation  breaks  down  (Westwater,  1972). 

A general  solution  to  the  problem  is  available  in  the  form  of  the  Mie  (1908) 
theory  solutions.  However,  the  extinction  efficiency  factor  must  be  computed 
for  each  drop  size  and  temperature  pair  at  a given  wavelength  and  then  inte- 
grated across  the  distribution  of  droplet  sizes  within  the  ensemble  in 
order  to  yield  the  extinction  coefficient.  This  calculation  would  involve 
significant  expense  if  implemented  within  the  intricacies  of  an  operational 
computation  scheme.  Thus,  there  is  a distinct  advantage  to  employing  an 
approximat ion . 

This  section  describes  a model  of  the  interaction  between  the  radia- 
tion field  and  the  precipitation  droplets  which  approximates  the  extinction 
due  to  the  presence  of  both  small  and  large  droplets  and  which  does  not 
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require  undue  amounts  of  computation  time.  It  should  be  emphasized  that  no 
attempt  was  made  to  develop  an  exact  radiative  transfer  model  in  which  the 
scattering  is  treated  as  a source  function.  Instead,  the  model  treats 
scattering  as  an  enhancement  in  the  extinction.  The  line  of  reasoning 
employed  has  been  to  extend  the  present  treatment  of  nonprecipitating 
e loud s to  include  the  effects  of  large  droplets.  In  this  sense  it  is  an 
analytical  extension  of  the  Staelin  approximation.  The  results  described 
provide  wavelength,  temperature,  and  precipitation  parameter  dependencies 
of  the  extinction  coefficient  in  a relatively  simple  analytical  form.  The 
dependence  of  the  extinction  coefficient  on  rainfall  rate  is  shown  to  com- 
pare favorably  with  comparisons  based  on  the  full  Mie  Theory  and  radar 
measurements.  In  addition,  comparisons  with  the  Staelin  approximation  are 
also  presented. 

2.2  Theoretical 


The  unit  volume  extinction  coefficient  due  to  an  ensemble  of  droplets 
is  given  by  the  convolution  of  the  extinction  cross  section  for  a droplet 
of  given  radius  with  the  drop  size  distribution  function: 
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is  unit  volume  extinction  coefficient  (neper  m ) 

is  drop  size  distribution  function  (cm  '**  ym 

is  the  Mie  efficiency  factor  for  extinction 

is  the  complex  index  of  refraction  (dimensionless) 

is  wavelength  (cm) 

is  drop  radius  (urn) 

is  absolute  temperature  (°K) 
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The  efficiency  factor  (which  is  a function  of  temperature,  droplet  radius, 
and  wavelength)  is  obtained  from  the  solution  to  the  problem  of  diffraction 
of  a plane  monochromatic  wave  by  a homogeneous  dielectric  sphere  of  arbi- 
trary radius  (Mie,  1908): 
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— Z [21  * 1)  Re  a.  (fl(X,T),q)  + b (fl(X,T) ,q) 

a 1=1  \ 


where  the  coefficients  and  in  the  Mie  series  are  obtained  by  the 

boundary  conditions  at  the  surface  of  the  sphere  and  at  infinity,  and  where 
4 

q = [2iTr/(X  x 10  )]  is  called  the  dimensionless  drop  size  parameter.  It  is 
a measure  of  the  phase  shift  across  the  droplet. 

The  complex  index  of  refraction  of  water  at  centimeter  wavelengths  is 
given  by  the  Debye  Formula  (1929)  : 

~ e - e 

~2  o °°  _ „ 

" = t --fryo  * ■ 2-> 

Here  c and  e are  the  static  and  optical  dielectric  constants  and  X is 
o 00  1 o 

the  Debye  relaxation  wavelength.  The  constants  are  temperature  dependent 
and  their  values  are  calculated  using  the  values  of  Grant  et  al  (1957)  . 

Since  the  complex  index  of  refraction,  n,  is  a function  of  temperature 
and  wavelength,  QE(n,X,r)  must  be  specified  for  each  X and  T,  as  well  as  for 
each  r.  Due  to  the  series  summation  in  Expression  2-2,  this  is  not  an  easy 
matter  computationally. 

2.3  Rayleigh  Limit  - Staelin  Approximation 


An  alternative  form  for  the  Mie  extinction  efficiency  factor  is  given 
by  a power  series  expansion  in  the  dimensionless  drop  size  parameter. 

QE(X,r,T)  = qZ (X ,T)  + q3  A(X,T)  + q4  B(X,T)  + ...  2-4 


where 


4 Im  {-K}  = 4 Im 


h2  - 1 
fi2  + 2 


~2  ~4  -2 

4 T /n  - 1 n + 27  n + 38 

- ? Im  -T—  — 72 — : 

\ n + 2 2n  + 3 


8/3  Re  ( n--~— -Y 

U * + 2) 
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This  expansion  is  valid  for  all  q,  but  is  especially  useful  in  the  limit 
where  q <<  1;  i.e.,  when  the  wavelength  of  interest  is  much  larger  than  the 
droplet  size.  At  microwave  wavelengths  this  holds  true  for  nonprecipitating 
c louds . 

In  this  limit,  the  Mie  solution  has  the  following  asymptotic  behavior: 


q+0;  QE-qZ(X,T). 


Rayleigh  Limit  (small  droplets) 


0 a 

For  nonprecipitating  clouds  the  average  droplet  radii  range  from  10  to  10“  pm. 

_3  -2 

At  a wavelength  of  1 cm,  the  q parameter  is  thus  ot  the  order  of  10  ' to  10  ; 

i.e.,  q <<  1 . This  falls  within  the  Rayleigh  limit  region. 

Substituting  into  Expression  2-1, 


Y 


L 


/ 


N (r)  qZ  (X  ,T)  Ttr2  dr 


mF(X,T) 
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where 


m 

and 


J N(r)  r^  dr 
o 


is  the  mass  density  of  the  droplet  ensemble 


F(X,T)  is  a function  of  wavelength  and  temperature  only.  Tins  function 
has  been  evaluated  empirically  by  Staelin  (1966)  for  which  the  expression 
becomes : 

. 0.0122  (291 -T) 

1 . 0016  x 10  m 10  (neper  m ) . 2-6 


This  says  that  the  unit  volume  extinction  of  a nonprecipitating  cloud 
is  proportional  to  its  mass  density  regardless  of  the  nature  of  the  droplet 
distribution.  This  is  not  true,  however,  for  a precipitating  cloud  and 
this  approach  is  not  applicable. 

2.4  The  ERT  Precipitation  Extinction  Approximation 

The  approximation  to  be  used  for  precipitating  clouds  is  based  on  the 
character  of  the  variation  of  Q (X,r,T)  with  r and  the  nature  of  the 
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relevant  drop  size  distribution  function  N(r).  The  upper  portion  of  F gure 
2-1  portrays  the  range  of  droplet  radii  for  fair  weather  cumulus,  cumulus 
congestus,  and  rain  with  a rainfall  rate  of  25.4  mm  hr  \ The  broken  line 
indicates  the  variation  of  (A,r,T)  with  r at  \ = 1.0  cm  and  T = 298°K. 
Note  that  while  the  size  distributions  of  nonprecipitating  clouds  fCu)  are 
limited  to  the  Rayleigh  region,  the  rainfall  distribution  extends  through 
the  transition  region  to  the  large  sphere  region  where  the  asymptotic 
limit  of  the  Mie  series  (Expression  2-2)  approaches  2: 


q 00  ; QL;  + 2. 


Large  Sphere  Limit 


An  average  raindrop  has  a characteristic  radius  on  the  order  of  10"^  pm. 

At  a wavelength  of  1 cm,  the  q parameter  is  of  order  1.  The  Rayleigh  limit 
is  no  longer  valid  and,  furthermore,  it  must  be  recognized  that  the  extinc- 
tion factor  will  vary  orders  of  magnitude  within  the  size  range  of  the  rain- 
fall droplets.  For  complete  accuracy,  the  full  Mie  Theory  using  Eq . 2-2 
should  be  used  in  this  range  of  drop  sizes.  This  would  be  feasible  if  all 
raindrops  had  the  same  radius  and  atmospheric  temperature  structure  was 
isothermal.  However,  distributions  of  droplet  radius  and  atmospheric  vertical 
temperature  structure  make  a layered  computation  using  Eq.  2-2  computationally 
complex . 

The  approach  taken  here  was  to  assume  that  the  Mie  expression  (Eq . 2-4) 
could  be  approximated  by  a small  number  of  terms.  Also,  the  number  of  terms 
included  was  assumed  to  be  dependent  on  the  size  range  of  the  droplets 
such  that  the  attenuation  coefficient  can  be  expressed  as: 


/ 


N(r)  QP  (r,A,T)  "r  dr 


/' 


N(r)  (r^)  irr  dr  + 


/2' 


(r)  Qe  (r2)  nr  dr 


N(r)  Qp  (r  ) irr2  dr 

fc  oo 
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n 

for  n + 1 size  intervals.  Ideally,  each  of  the  above  integrals  should  be 
expressible  in  closed  form. 
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Based  on  QR  (X,r,T)  vs.  r plots  such  as  in  Figure  2-1,  three  representative 
regions  were  selected:  ’ 


(a)  Region  I 

(b)  Region  II 

(c)  Region  III 


Rayleigh  limit  region 
Transition  region 
Large  sphere  region 


< r„ 


rR  1 r 1 rc 

r < r 
c 


The  approximations  to  the  Mie  efficiency  factors  in  the  three  regions  are  as 
fo 1 lows : 


QE 


qZ 

qZ  + Aq^  + Bq4 


for  0 < r <_  rR 

for  i-  < r < r 
R — — c 

for  r > r 
— c 
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The  quantities  q,  Z,  A, and  B are  as  previously  defined. 

The  limits  r^  and  rR  are  specified  in  the  following  manner.  The  upper 
limit  to  the  Rayleigh  region  is  taken  where  the  dimensionless  drop  size 

7 

parameter  reaches  0.10.  Tins  is  approximately  the  point  at  which  the  q' 
term  in  iiq.  2-4  reaches  10%  of  the  leading  term  in  q.  Thus, 


qR  = 0.10 

and 

rR  = (1.59  x 10  2)X. 

The  lower  limit  of  the  large  sphere  region  is  located  by  specifying 
the  position  of  the  optical  resonance  region.  The  general  criteria  adopted 
is : 


4 

X x 10 


where 

| n | > r 

The  fit  of  this  approximation  to  the  Mie  Theory  is  demonstrated  in  Figure 
2-1  and  is  quite  reasonable. 
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Within  the  transition  region  bounded  by  (r^>rc),  the  drop  size  parameter 

takes  on  values  in  the  range  of  0.1  < q < 1.0.  The  higher  order  terms 
3 4 

Aq*  and  Bq  , therefore,  contribute  significantly  to  the  magnitude  of  , 
while  higher  orders  in  q are  small  corrections. 

Based  on  this  development,  Expression  2-7  is  truncated  to  a series  of 
three  integrals: 

/R 

Y - J qZN(r)  nr 
o 


dr  + 


i 

/ 


(Aq' 


4 

Bq  ) 


N(r)  nr  dr 


+ 


2 N (r ) 


nr 


dr 
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where  all  the  terms  are  as  previously  defined. 


2.5  Drop  Size  Distribution  Function 

In  order  to  evaluate  Eq.  2-9  in  closed  form,  the  drop  size  distribution 
must  be  specified.  The  particular  form  of  the  drop  size  distribution  N(r) 
is  empirically  determined  and  for  the  purpose  of  our  study,  the  most  general 
specification  for  cloud  droplet  spectra  is  the  Deirmendjian  (1964)  distribu- 
tion : 

C C 

N(r)  = Ar  * exp  {-Br  “}  2-10 


which  characterizes  the  ensemble  in  terms  of  four  parameters,  two  of  which 
(A  and  B)  are  scale  parameters,  the  others  being  shape  parameters.  Fair 
weather  and  stratus  clouds  conform  nicely  to  this  characterization.  A 
simple  form  of  the  Deirmendjian  distribution  may  be  used  to  describe  p*ecipi- 
tat ion . 

A review  of  the  literature,  especially  that  pertaining  to  radar  studies 
of  precipitation,  indicated  that  the  Marshal  1 -Palmer  (1948)  model  of  drop 
size  distribucion  for  precipitation  is  the  most  appropriate  for  studies  of 
satellite-measured  brightness  temperatures.  The  virtues  of  this  model  are 
that:  (1)  an  empirical  expression  exists  which  relates  the  drop  size 
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distribution  (useful  in  computations  of  radiati'  ■;  interaction)  to  the  instan- 
taneous rainfall  rate  (useful  in  synoptic  studies),  and  (2)  the  form  of  the 
expression  is  essentially  that  of  the  Deirmendjian  expression  used  in  modeling 
cloud  droplet  distributions.  The  Marsha  1 1 -Pa lmer  model  is  expressed  as: 

"^3 

N(r,R)  = exp  - (k0  rR  } 2-11 

where 

k = 8.0  x 10'6 

k2  = 8.156  x 10-3 

k3  = 0.21 

and 

N(r,R)  is  the  number  density  of  particles  per  unit  size  range  (cm  ''pm*) 

r is  the  droplet  radius  (pm)  for  a given  rain  rate. 

R is  the  rainfall  rate  (mm/hr). 

This  is  essentially  the  Deirmendjian  distribution  with  parameters: 


B = k_,R  J 

fL 

s - 0 

C;  = 1 

Substituting  bq.  2-11  into  bq . 2-9  and  evaluating,  the  extinction  coefficient 
becomes  a function  of  wavelength,  temperature,  and  rainfall  rate: 

y ( A , T , R)  = E C F (A)  G ( A , T ) H ( A , R)  R°’21n  2-12 

. nn  n n 

n=l 

where 

C,,  are  constants 
N 

and 

3re  the  indicated  variables. 
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J.b  Discussion  of  the  Approximation 


The  derived  expression  for  the  extinction  coefficient  y(,\,T,R)  is 

essentially  a power  series  in  the  rainfall  rate,  R,  with  lowest  term  going 

0 "■  1 1 47 

as  R ' “ and  highest  term  going  as  R ' . Temperature  and  wavelength 

dependence  are  given  by  the  coefficients  in  this  finite  scries.  It  can 
be  shown  that  this  solution  reduces  to  the  Staelin  approximation  in  the  limit 
as  the  rainfall  rate  approaches  cero.  This  limit  is  demonstrated  in  Figure 
2-2  giving  the  dependence  of  extinction  coefficient  on  temperature  at  low 
rainfall  rates  (R  = 0 and  1 mm  hr  and  at  an  intermediate  rainfall  rate 
(R  = 39.5  mm  hr  *).  Note  that  althcugh  the  intermediate  rates  are  essenti- 
ally temperature  independent,  there  is  a marked  variation  (by  a factor  of 
4 to  5)  at  lower  rainfall  rates.  Additionally,  the  extinction  due  to  the 
cloud  alone  (R  = 0)  is  typically  an  order  of  magnitude  less  than  that  due 
to  medium  precipitation  (R  = 39.5  mm  hr  ^).  The  importance  of  the  cloud 
contribution  decreases  with  increasing  temperature.  Comparisons  between  the 
approximation  and  both  Mie  Theory  and  the  Staelin  Formula  (Rayleigh  scatter- 
ing) can  be  made.  Various  authors  have  made  full  Mie  computations  at  radar 
wavelengths  based  on  convolution  with  the  Marsha  1 1 -Pa lmer  distribution. 

Wexler  and  Atlas  (1963)  and  Shifrin  and  Chernyak  (1968)  have  demonstrated 
that  the  volume  extinction  coefficient  is  a near-linear  function  of  rainfall 
rate.  This  behavior  is  reproduced  by  the  approximation.  Figure  2-3  demon- 
strates the  extinction  vs.  rainfall  rate  relationship  at  i = 1.35  cm,  T - 283°K 
for  full  Mie  Theory,  the  precipitation  approximation,  and  the  Staelin  Formula. 
The  extreme  underestimation  of  the  extinction  by  the  Staelin  approximation 
is  due  to  the  total  neglect  of  scattering  and  higher  order  absorption  by 
the  Rayleigh  approximation.  The  prec ipitat  ion  approximation  includes  these 
effects  to  first  order  and  improves  to  within  85%  of  the  full  Mie  series 
solution  for  the  worse  case  at  60  mm/hr.  Figure  2-4  compares  three  radar 
extinction  vs.  rainfall  rate  results  with  the  approximation.  On  the  basis 
of  these  calculations  it  is  concluded  that  the  approximation  significantly 
improves  the  rainfall  rate  dependence  of  the  extinction  coefficient  predicted 
by  the  Staelin  Formula,  approaching  the  behavior  of  the  full  Mie  series 
solut  ion . 

The  frequency  dependence  of  the  extinction  coefficient  is  significant 
in  determining  information  retrievable  from  mult i spectral  analysis.  Figure 


13 


2-5  demonstrates  the  frequency  dependence  of  the  ERT  approximation  as 
compared  to  the  Staelin  Formula  (Rayleigh  Scattering -only  curve)  and  two 
full  Mie  Theory  calculations:  (1)  a Deirmend j ian  distribution  (parameters 

-3 

C = 5.0,  = 0.6)  with  a mode  radius  of  400  um  and  a mass  density  of  1.0  g m ; 

*■  _ 3 

and  (2)  a Marshal  1 -Pa lmer  distribution  with  a mass  density  of  1 g m ' . Figure 
2-5  demonstrates  that  the  approximation  models  the  frequency  dependence  of 
full  Mie  Theory  calculations  quite  well. 
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Extinction  Coefficient  (Neper/cm) 


Mass  radius  -400  a ~ 

Computed  mass  density  - / 00 a g m~ 5 

Cj-5  0 Ideirmendjian  distribution 
A-  0^0.6/ precise  computation 
— ERT  Approximation 
^-Marshall -Palmer  distribution 
lg  m~3  ( precise  computation) 

— Rayleigh  Scattering  only 


Frequency  (GHz) 


Figure  2-5  Frequency  Dependence  of  the  Extinction  Coefficient  for 
Various  Models 


3.  TEMPERATURE  PROF ILINC  ANALYSIS 


3 . 1 Scope 

Aspects  of  atmospheric  temperature  profiling  by  passive  microwave 
techniques  were  systematically  analyzed.  These  included  the  following: 

• Hffeets  of  line  shape  model  on  the  computation  of  weighting 
functions  and,  therefore,  simulated  brightness  temperatures. 

• Comparison  of  the  use  of  the  118  GHz  0o  complex  in  tempera- 
ture profile  with  that  of  the  60  GHz  0^  complex. 

e Evaluation  of  the  effects  of  clouds  and  precipitation  on 
temperature  profiling  in  the  two  spectra]  regions. 

The  results  of  these  analyses  are  presented  in  this  section. 

3.2  Approacti 

3.2.1  Line  Shape  Analysis 

Previously,  the  ERT  simulation  programs  made  use  of  the  Van  Vleck- 
Weisskopf  line  shape  model.  For  this  study,  the  Gross-Reber  model  was 
adopted.  Evaluations  of  this  change  on  temperature  profiling  simulation 
results  were  performed  by  comparing  the  differences  in  absorption  computed 
using  the  two  models.  Computations  were  made  at  a number  of  frequencies  and 
at  a number  of  pressures  and  concentrations.  In  addition,  typical  weighting 
functions  were  computed  using  the  models  and  compared. 

3.2.2  Temperature  Profiling  Analysis 

The  remaining  two  tasks  identified  previously  were  accomplished  by 
simulation  inversion  analysis  using  ERT  programs.  The  inversion  procedure, 
the  statistical  method  for  estimation  of  parameters,  has  been  well  documented 
in  previous  reports  such  as  Gaut , Reifenstein  and  Chang  (1972).  This  proce- 
dure is  shown  schematically  in  Figure  3-1.  Briefly,  the  approach  uses 
physical  and  interaction  models  operating  on  an  ensemble  of  atmospheric 
data  (including  clouds  and  precipitation)  to  create  a correlation  matrix 
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Preceding  page  blank 
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Error  ino/ysis 

I 

Figure  3-1  Flow  Diagram  for  Application  of  the  Data  Inversion  Methodology  to 
Parameter  Estimation  for  a Radiometric  Experiment.  (The  and  d_ 
data  sets  may  be  simulated  for  study  purposes,  as  shown  by 
dashed  lines.) 


or  "1>-Matrix"  between  brightness  temperatures  at  selected  frequencies  and 
geophysical  parameters  of  interest.  This  so-called  1)-Matrix  is  then  used 
to  "infer"  from  actual  radiometric  measurements  the  values  of  the  correspond- 
ing geophysical  parameters.  In  the  analysis  and  evaluation  made  where  real 
radiometric  measurements  are  not  available,  a simulation  procedure  is  used 
to  generate  a set  of  "experimental  data."  The  simulation  procedure  used  is 
identical  to  that  used  to  generate  the  data  for  t It e D-Matrix.  However,  a 
separate  and  different  ensemble  of  atmospheric  and  surface  data  is  used 
for  this  evaluation  data  set. 

for  the  purposes  of  this  study  actual  historic  radiosonde  data  were 
used  to  describe  the  vertically  correlated  temperature  structure  of  the 
atmosphere . Similarly,  radiosonde  data  were  used  to  specify  the  vertical 
distribution  of  atmospheric  water  vapor.  Since  radiosonde  data  do  not 
contain  basic  information  on  the  presence  and  types  of  clouds,  these  have 
to  be  inserted  according  to  some  realistic  procedure.  Statistics  of  the 
probability  of  occurrence  for  each  of  six  major  cloud  categories  (including 
clear-shy  conditions)  corresponding  to  the  season  and  locations  of  the 
radiosonde  data  were  used,  together  with  a random  number  generator,  to  insert 
the  proper  distribution  of  clouds  ir.co  the  ensemble  of  radiosonde  data, 
bach  of  the  major  cloud  types  lias  a physical  description  in  terms  of  its 
vertical  extent,  drop  sice  distribution  and  liquid  water  content.  These 
parameters  are,  in  fact,  the  ones  whicli  get  inserted  into  the  radiosonde 
profiles  to  create  the  ensemble  of  atmospheric  dat  i . 

In  some  of  the  analyses  performed  for  this  study,  a different  approach 
was  adopted  in  creating  ensembles  of  atmospheric  data  for  evaluation.  In 
order  to  clearly  identify  the  effects  of  clouds  and  precipitation  on  inver- 
sion results,  different  ensembles  of  atmospheric  data  were  created,  using 
the  same  set  of  temperature  profiles,  to  represent  conditions  with: 

• Clear  sky 

e A cloud  with  fixed  thickness  and  liquid  water  content 

• Precipitation  with  fixed  liquid  water  content 

Each  of  these  ensembles  were  "inverted"  using  the  same  D-Matrix  that  was 
created  using  cloud  and  precipitation  conditions  representative  of  the 
season  and  station  location.  By  comparing  the  inversion  results  either  in 
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terms  of  the  ensemble  statistical  rms  errors  or  in  terms  of  individual 
profile  accuracies,  the  effects  of  clouds  and  rain  (and,  therefore,  drop 
size)  can  be  isolated. 

3.3  The  Attenuation  Due  to  Oxygen 

The  attenuation  of  microwave  radiation  near  frequencies  of  60  GHz  and 
118  Cllz  is  due  primarily  to  the  resonance  of  oxygen  molecules.  These 
resonances,  caused  by  magnetic  dipole  rotational  transitions,  occur  at 
several  frequencies  within  the  50  to  60  GHz  range  and  at  118.7505  GHz. 

Each  resonance  line  has  a line  width  related  to  the  radiation  transition 
lifetime  which  is  broadened  by  molecular  collisions  in  a manner  which  is 
pressure  and  temperature  dependent.  As  a result,  the  individual  lines 
between  50  and  60  GHz  are  clearly  defined  at  very  low  pressures,  but  combine 
to  form  a continuous  absorption  region  at  pressures  of  one  atmosphere. 

To  compute  the  microwave  absorption  coefficient  due  to  oxygen,  it  is 
necessary  to  sum  the  absorption  caused  by  each  rotational  state  found  in 
appreciable  populations  at  atmospheric  temperatures.  This  requires  Know- 
ledge of  the  line  shapes  and  line  width  parameters,  information  which  is 
currently  provided  from  theory  but  must  be  empirically  supported.  Previous 
studies  at  ERT  have  followed  the  approach  presented  by  Meeks  and  Li 1 ley 
(1963)  using  their  line-width  expression  with  the  line-shape  factor  defined 
by  Van  Vleck  and  Weisskopf  (1945)  to  determine  the  contribution  of  each 
resonance  line  to  the  absorption  coefficient.  However,  recent  research  (e.g., 
Reber,  1972),  has  indicated  that  the  form  factor  developed  by  Gross  (1955) 
agrees  more  accurately  with  measurements  of  microwave  attenuation.  Hie  Gross 
line  shape  was  thus  adopted  for  this  study  and  comparisons  were  made  between 
the  attenuation  predicted  by  the  Van  Vleck-Weisskopf  expression  and  by  the 
Gross  form  factor. 

3.3.1  The  Oxygen  Line  Shapes 

The  theory  behind  collision  broadening  of  spectral  lines  was  first 
developed  by  Lorentz  (1909)  for  resonant  phenomena  at  optical  wavelengths 
and  by  Debye  (1929)  for  nonresonant  phenomena.  Van  Vleck  and  Weisskopf 
(1945)  analyzed  and  reformulated  their  theories  to  derive  a general  formula 
for  the  absorption  coefficient  based  upon  an  assumption  of  strong  molecular 


collisions.  As  a result  of  their  research,  the  following  line-shape  expres- 
sion was  derived: 


( v,  v ) = 
n 


where 


Av  __  Av 

2 T + r~=i  ~ 

(v  - v)“  + Av“  (v  + v)~  + Av" 
n n 
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v = observation  frequency  (GHz) 

v = resonant  frequency  of  the  nth  rotational  transition  (GHz) 

Av  = line  width 


Gross  (1955)  reexamined  all  of  the  above  theories  and  the  assumptions 
upon  which  they  were  based.  In  the  Van  V 1 eck-Wei sskopf  approach,  molecular 
collisions  were  assumed  to  leave  the  oscillator  positions  and  velocities 
distributed  according  to  the  Boltzmann  distribution  appropriate  to  the 
applied  field  at  the  time  of  collision.  Gross  showed  that  this  required 
unrealistically  large  oscillator  velocities  and  that  the  theory  failed 
at  zero  and  infinite  frequencies,  lie  thus  modified  the  theory  by  assuming 
that  the  oscillator  positions  are  unchanged  as  a result  of  collision,  but 
that  the  velocities  are  altered  to  the  Boltzmann  distribution  corresponding 

to  the  time  of  impact.  This  led  to  a reformulation  of  the  form  factor: 

*> 

4v  Av 

^ (v»vn)  “ 2 2 ~ 2 2” 

(v  “ - v")“  + 4v"Av" 
n 

where  the  variables  are  defined  above. 

Rebcr  (1972)  made  an  extensive  series  of  measurements  of  zenith  atmo- 
spheric attenuation  in  the  50  to  70  GHz  range  and  compared  his  results  with 
those  predicted  by  the  Lorentz,  Van  Vleck-Keisskopf  and  Gross  form  factors, 
fitting  the  theoretical  values  to  his  data,  he  found  that  the  Gross  and 
Lorentz  line  shape  resulted  in  an  rms  error  of  8.24%  and  8.25%,  respectively. 
The  Van  V leek- Wei sskopf  line  shape,  cn  the  other  hand,  showed  an  rms  error 
of  15.70%.  He  also  showed  that  the  Van  Vleck-Weisskopf  form  factor,  while 
generally  agreeing  with  the  others  between  50  and  60  GHz,  skews  the  total 
attenuation  curve,  especially  above  60  GHz,  and  that  at  70.26  GHz,  it  pre- 
dicts 30  to  40%  more  zenith  attenuation  than  his  measured  values  showed. 


3.3.2  A Comparison  Between  the  Van  Vleck-Weisskopf  and  the  Cross 
Line  Shapes 

To  determine  the  effect  of  the  differences  between  the  Van  V lech -Wei ss- 
kopf  and  Cross  line  shapes  on  the  simulation  of  atmospheric  attenuation, 
computations  of  the  oxygen  absorption  coefficient,  y,  were  made  at  1 GHz 
intervals  from  30  to  130  GHz  for  pressures  representative  of  the  atmosphere. 
These  arc  shown  in  Figure  3-2  for  the  Van  Vleck-Weisskopf  line  shape  and  in 
Figure  3-3  for  the  Cross  line  shape.  Excellent  agreement  is  seen  near  b0  til Iz 
and  119  GHz  where  the  resonant  frequencies  occur,  since  the  two  formulations 
agree  in  value  on  the  oxygen  lines.  It  thus  appears  that  for  near-resonance 
microwave  channels,  either  line  shape  provides  a good  absorption  prediction. 

Away  from  the  oxygen  lines,  significant  differences  are  found  between 
the  two  expressions.  Below  35  GHz  the  Van  Vleck-Weisskopf  expression  esti- 
mates attenuation  above  that  predicted  by  the  Gross  expression,  while  between 
35  GHz  and  50  GHz  it  predicts  coefficients  lower  than  those  predicted  by 
Gross,  although  remaining  within  10%  of  the  value.  More  pronounced  differ- 
ences are  found  between  65  GHz  and  110  GHz,  since  in  this  frequency  range 
the  Van  Vleck-Weisskopf  absorption  coefficient  far  exceeds  the  Gross  one. 

This  is  especially  true  above  85  GHz  where  the  coefficients  disagree  by  a 
factor  of  2.  These  results  are  in  line  with  the  findings  of  Kebcr  (1972] 
discussed  above  and  indicate  that  the  use  of  the  Gross  line  shape  for  non- 
resonant  frequencies  is  of  great  importance  to  an  accurate  simulation  of 
atmospheric  attenuation. 

Figures  3-4  and  3-5  show  the  weighting  functions  computed  from  the  two 
expressions  for  frequencies  of  55.29  GHz  and  54.3  GHz,  respectively . These 
were  generated  by  the  simulation  program  for  a tropical  atmosphere  over  a 
calm  ocean.  Both  of  these  frequencies  arc  near  line  centers  and,  as  would 
be  expected  from  the  analysis  given  above,  show  no  significant  differences. 
Similar  agreement  was  found  for  other  channels  in  the  50  to  60  GHz  range 
and  near  118  GHz,  indicating  that  the  inversion  for  atmospheric  temperature 
should  provide  the  same  results  using  cither  the  Van  Vleck-Weisskopf  or  the 
Gross  line  shape.  However,  temperature  inversions  are  often  performed 
simultaneously  with  inversions  for  water  vapor  or  liquid  water  using  channels 
away  from  the  oxygen  lines.  Thus,  the  Gross  line  shape  was  incorporated 
into  the  INVERT  programs  to  provide  t he  maximum  accuracy  in  the  simulation 
oi  radiative  transfer. 
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3.4  Inversion  Analysis  Results 


3.4.1  Selection  of  Spectral  Channels 

The  starting  point  of  the  analysis  was  the  set  of  five  channels  in  the 
60  GHz  complex  previously  selected  by  F:RT  as  appropriate  for  temperature 
profiling  in  the  troposphere  and  lower  stratosphere.  These  channels  are 
shown  in  Figure  3-6  in  terms  of  their  weighting  function  computed  for  mid- 
latitude  (Figure  3-6a)  and  for  tropical  (Figure  3-6b)  clear-sky  conditions 
with  an  ocean  background.  The  peaks  of  the  weighting  functions  are  locate’ 
at  approximately  2,  7,  13,  17  and  23  km  and  do  not  appear  to  differ  signifi- 
cantly between  the  tropical  and  mid  latitude  conditions.  The  shape  of  the 
weighting  function  for  the  58.80  GHz  channel  with  peak  value  at  17  km  appears 
compressed  with  a higher  maximum  in  the  tropical  case  due  to  the  effect  of 
the  tropical  tropopause  which  is  located  at  approximately  this  altitude. 

Comparable  channels  were  sought  in  the  wing  of  the  118  GHz  line,  where 
difficulty  was  encountered  in  defining  a channel  with  a peak  value  corre- 
sponding to  that  of  the  53.29  GHz  line;  i.e.,  with  peak  at  about  2 km.  The 
selected  set  of  frequencies  and  their  corresponding  weighting  functions  is 
graphically  illustrated  in  Figure  3-7,  with  the  midlatitude  ocean  case  shown 
in  Figure  3-7a  and  the  tropical  ocean  case  in  Figure  3-7b.  Compared  to  the 
corresponding  weighting  functions  for  the  60  GHz  channels,  it  is  clear  that 
with  tl'.e  exception  of  the  117.76  GHz  channel,  the  match  between  the  two  sets 
of  weighting  functions  is  quite  good.  Further  analysis  of  the  117.76  GHz 
channel  suggests  that  the  difficulty  in  defining  a channel  with  peak  at 
low  altitudes  is  due  to  the  strong  effects  of  water  vapor  attenuation  con- 
tributed by  the  183  GHz  water  vapor  line.  Figure  3-8  shows  the  weighting 
functions  for  the  same  set  of  frequencies,  but  computed  for  a dry  atmo- 
sphere;  i.e.,  with  integrated  water  vapor  of  0 g cm*’.  Under  this  condition, 
the  117.76  GHz  channel  peaks  at  approximately  2 km  and  the  change  in  the 
118  GHz  channel  is  also  significant.  Changes  in  the  weighting  functions 
for  the  other  channels  are  less  apparent  due  to  the  fact  that  the  water 
vapor  in  the  atmosphere  occurs  mainly  at  levels  below  5 km. 
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Temperature  Inversion  Results:  The  General  Case 


A set  of  general  inversion  analyses  was  performed  using  the  60  GHz  and 
the  118  GHz  channels  separately,  with  cloud  statistics  appropriate  to  the 
selected  regions  and  seasons.  Also  included  in  these  analyses  were  channels 
selected  to  probe  for  water  vapor  and  cloud  liquid  water.  Table  3-1  sum- 
marizes the  two  sets  of  frequencies  used.  In  all  of  the  downlooking  cases, 
the  channels  were  nadir  viewing.  In  the  uplooking  cases,  the  instruments 
sensed  directly  at  zenith. 


TAdLL  3-1 

FREQUENCIES  SELECTED  FOR  ANALYSIS  OF  THE  GENERAL  CASE 


Purpose 

v < 60  (GHz) 

v > 117  (GHz) 

Temperature  Profiling 

53.29 

117.76 

54.30 

118.00 

55.46 

118.22 

58.80 

118.5 

59.50 

118. 675 

Water  Vapor  and  Liquid  Water 

18.0 

18.0 

21.0 

21.0 

22.235 

22.235 

31.6 

31  .6 

Inversion  analysis  was  performed  for  two  climatic  conditions:  midlati- 
tude spring  and  tropical  summer.  Spring  in  the  midlatitudes  is  a season 
of  rapid  and  frequent  changes  in  air  masses,  and  thus  in  atmospheric  temper- 
ature and  moisture  structures.  Data  from  this  climatic  regime  is  thus 
representative  of  a high  degree  of  atmospheric  variability.  Climatic  condi  • 
tions  in  the  tropics,  on  the  other  hand,  are  quite  stable  and  the  natural 
variability  in  both  temperature  and  moisture  is  quite  small. 

The  inversion  results  for  the  downward- looking  cases  are  shown  in 
Figures  3-9  to  3-12,  expressed  as  rms  errors.  The  midlatitude  analysis 
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was  performed  using  a D-Matrix  generated  from  84  atmospheric  profiles  and 
evaluated  against  36  simulated  measurement  sets.  The  tropical  analysis  used 
a D-Matrix  generated  from  76  atmospheric  profiles  and  evaluated  against  36 
simulated  measurements  sets.  In  each  case,  the  data  plotted  represent  the 
rms  difference  for  the  36  evaluation  samples  between  the  inversion  results 
and  the  input  profile  values.  Also  shown  are  the  expected  a priori  vari- 
ability  of  each  of  the  parameters  at  each  level  in  the  atmosphere.  These 
a priori  values,  computed  from  the  data  used  to  create  the  D-Matrices,  can 
be  interpreted  as  the  rms  error  which  would  result  if  the  profiles  were  simply 
"guessed  at"  based  on  climatology  alone.  The  differences  between  the  a priori 
variability  and  the  inversion  rms  error  is  a measure  of  the  improvement  of 
the  results  inferred  from  radiometric  measurements  over  those  estimated 
from  climatological  data. 

In  examining  the  figures,  the  differences  in  the  variability  between 
the  tropical  and  midlatitude  spring  atmospheres  are  clearly  reflected  in 
the  a priori  statistics.  In  the  tropospheric  levels  the  rms  variability 
of  temperature  for  the  midlatitude  case  is  approximately  7°K,  while  for  the 
same  altitude  region,  the  rms  variability  of  the  tropical  atmosphere  is 
only  of  the  order  of  1°K.  A similar  difference  exists  for  the  moisture 
profiles,  especially  at  the  lower  levels.  Because  of  this  difference  in 
a priori  variability,  it  can  be  expected  that  remote  sensing  data  are  more 
useful  in  the  midlatitude  zones.  The  results  of  the  error  analysis  shown  in 
these  figures,  in  fact,  demonstrate  this  phenomenon.  For  the  midlatitude 
situation,  either  of  the  channel  combinations  reduces  the  rms  error  in 
estimating  temperature  profiles  by  some  5°K  to  6°K  over  climatology  in  the 
tropospheric  layer  above  2 km,  while  the  differences  between  the  inversion 
rms  error  and  a_  priori  statistics  for  the  tropical  temperature  case  is 
only  near  1°K. 

In  comparing  the  midlatitude  results  obtained  using  the  higher  fre- 
quency oxygen  channels  with  those  using  the  60  GHz  channels,  it  is  evident 
that  at  levels  below  6 km  the  higher  frequency  channels  are  not  as  effec- 
tive in  measuring  temperature.  This  difference  is  probably  a direct  result 
of  the  fact  that  the  weighting  functions  for  the  set  of  channels  in  the 
118  GHz  line  were  shown  to  be  sensitive  to  atmospheric  water  vapor. 
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Above  0 km,  whore  water  vapor  exists  at  highly  reduced  concentration,  no 
substantial  difference  between  the  two  sets  of  channels  is  evident. 

Inversion  error  analysis  was  also  performed  for  the  upward- looking 
case  using  the  same  sets  of  channels.  Computations  were  made  only  for  the 
midlatitude  spring  climatic  condition.  The  results  are  shown  in  figure  3-13 
and  figure  3-14.  Since  these  channels  were  not  in  any  sense  optimized  for 
the  upward-looking  situations,  it  is  surprising  to  note  the  relatively  good 
performance  of  the  60  GHz  channels  in  the  lower  and  midtroposphere  with 
inversion  rms  errors  in  temperature  of  < 1.5°K  up  to  6 km.  The  relatively 
poorer  performance  of  the  118  GHz  channels  in  tropospheric  temperature 
probing  is,  again,  evident  in  the  upward- looking  case. 

3.4.3  effects  of  Clouds  and  Kain  on  Temperature  Profiling 

The  results  discussed  in  Se  t ion  3.4.2  are  based  on  a D-Matrix  with 
atmospheric  profiles  which  inelvde  clouds  where  the  frequency  of  occurrence 
corresponds  to  the  climatological  conditions  appropriate  to  the  location 
and  season.  Similarly,  the  evaluation  data  sets  of  36  profiles  also  con- 
tained cases  with  various  cloud  types.  However,  the  evaluation  sample  was 
too  small  to  permit  a systematic  evaluation  of  the  effects  of  clouds  on 
inversion  result  accuracies,  so  new  data  sets  were  created  with  clouds  and 
precipitation  of  known  properties  systematically  inserted  into  the  tempera- 
ture profiles.  This  made  it  possible  to  isolate  the  effects  of  cloud  thick- 
ness, cloud  liquid  water  content,  and  cloud  height  on  temperature  incursion 
accurac i cs . 

As  a preliminary  to  the  actual  inversion  analysis,  computations  were 

made  of  t lie  weighting  functions  of  selected  channels  in  the  presence  of 

hydrometeors  for  the  tropical  atmosphere  case.  Figures  3-15  and  3-16  show 

the  weighting  functions  for  the  two  lower  frequency  channels,  for  both 

sets  of  oxygen  channels,  in  the  presence  of  a nonprecipitating  cloud  of 

the  same  cloud  liquid  water  density  but  at  two  different  cloud  depths.  In 

_3 

the  first  instance  (Figure  3-15),  a uniform  cloud  of  0.6  g m ' density  was 
inserted  into  the  atmosphere  between  the  surface  and  4 km . In  the  second 
instance  (Figure  3-16),  a similar  cloud  having  the  same  density  and  base 
was  used  but  the  cloud  extended  to  8 km. 


a / 
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The  presence  of  these  clouds,  even  in  the  case  of  the  4 km  cloud,  is 
clearly  reflected  in  the  117.76  and  118.0  channels  where  the  attenuation 
due  to  liquid  water  and  the  associated  water  vapor  limits  cloud  penetration 
to  one  kilometer  Much  less  of  an  effect  is  seen  in  the  lower  frequencies, 
although  the  8 km  cloud  causes  substantial  changes  in  the  weighting  functions 
for  the  53.29  and  54.3  GHz  channels,  respectively.  For  the  4 km  cloud  the 
effect  on  the  53.29  GHz  channel  is  clear;  however,  the  54.3  channel  does  not 
appear  to  be  greatly  influenced  by  this  cloud. 

The  attenuation  due  to  an  ice  layer  above  a cloud  layer  is  shown  in 
Figures  3-17  and  3-18  where  an  ice  layer  2 km  thick  was  added  to  the  4 km 
cloud,  and  a 4 km  thick  ice  layer  was  added  to  the  8 km  cloud.  In  comparing 
these  figures  with  Figures  3-15  and  3-16,  it  is  clear  that  the  effect  of  ice 
on  the  oxygen  channels  is  minimal,  even  for  a 4 km  ice  layer. 

The  effects  of  precipitation  on  the  weighting  functions  are  shown  in 

Figures  3-19  and  3-20.  In  these  analyses  the  4 km  cloud  has  been  replaced 

by  a 4 km  layer  of  uniform  rain,  and  the  8 km  cloud  by  8 km  of  uniform 

rain  capped  by  a 4 km  layer  of  ice.  The  first  case,  that  of  the  4 km  rain 

layer,  is  representative  of  tropical  conditions  with  precipitation  being 

initiated  at  the  freezing  level.  The  latter  case  is  not  meant  to  be  realistic 

but  is  used  simply  as  a limiting  case  for  analysis.  In  comparing  these  two 

figures  with  those  for  the  clouds  of  comparable  thickness,  the  increased 

attenuation  effects  of  precipitation  are  clearly  evident.  These  effects  are 

due  to  the  enhanced  effects  of  scattering  as  modeled  in  the  precipitation 

approximation  (see  Section  2)  and  not  to  a difference  in  liquid  water  content. 

In  fact,  the  liquid  water  content  for  the  4 km  cloud  exceeds  that  of  the 
2 2 

rain  (0.24  g/cm  for  the  4 km  cloud  vs.  0.12  g/cm  for  the  4 km  rain).  The 
enhanced  attenuation  is  much  more  dramatic  for  the  lower  frequency  pair  of 
channels  since  for  the  4 km  cloud  case  attenuation  at  53.29  GHz  was  not 
significant.  With  precipitation  the  surface  becomes  obscured  to  all  four 
channels. 

These  results  also  clearly  demonstrate  the  frequency  dependence  of 
the  differences  between  the  ERT  precipitation  model  and  the  Staelin  approxi- 
mation for  clouds.  As  shown  previously  in  Figure  2-5,  the  difference  be- 
tween the  two  approximations  is  near  maximum  at  50  GHz.  The  differences 
in  weighting  function  at  this  frequency  range  should  be  large  between 
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those  computed  with  cloud  and  the  one  computed  with  rain.  At  higher  frequen- 
cies, the  absorption  coefficient  computed  by  the  two  different  approximations 
conveiges.  Differences  in  the  weighting  function  for  clouds  and  for  rain 
should  therefore  converge  to  zero  at  the  higher  frequencies,  and  appear  to  be 
doing  so  in  the  figures  depicting  the  118  GHz  line. 

The  presence  of  an  ice  layer,  even  of  4 km  thickness,  does  not  appear  tc 
have  any  noticeable  effect  on  the  weighting  functions. 

Inversion  analyses  were  systematically  performed  to  evaluate,  quanti- 
tatively, the  effects  of  hydrometeors  on  temperature  profiling.  The  results 
of  these  analyses  are  shown  in  Figures  3-21  and  3-22  for  channels  in  the  6U 
GHz  set,  and  in  Figures  3-23  and  3-24  for  the  corresponding  channels  in  the 
118  GHz  set.  In  each  case  the  oxygen  channels  were  used  in  conjunction  with 
the  water  vapor  and  surface  channels  defined  in  Table  3-1.  In  every  case  the 
D-Matrix  used  for  the  inversion  contained  the  distribution  of  clouds  appropri- 
ate to  the  climatic  regime. 

Figure  3-21  shows  the  inversion  results  obtained  for  the  60  GHz  channels, 
using  an  evaluation  sample  consisting  of  36  clear-sky  temperature  profiles;  the 
same  temperature  profiles,  but  with  a uniform  cloud  of  2 km  thickness  and  a 
density  of  1 g/m  , and  the  same  temperature  profiles  but  with  a 2 km  layer  of 

3 

rain  having  a uniform  liquid  water  density  of  1 g/m  . The  results  are  as 
one  would  expect  from  the  analysis  of  the  weighting  functions.  In  the  pres- 
ence of  clouds,  the  inversion  errors  increase  over  those  for  clear  sky  con- 
ditions. These  errors  are  larger  still  for  the  rain  case.  In  the  atmospheric 
layer  below  2 km  the  increase  in  rms  error  is  over  twofold  in  the  presence 
of  rain.  The  unexpected  feature  of  these  results  is  the  fact  that  the  increases 
in  rms  error  are  reflected  at  levels  of  the  atmosphere  significantly  above 
the  top  of  the  hydrometeor  layer.  At  4 km,  which  is  2 km  above  the  hydro- 
meteor layer,  the  rms  error  for  the  precipitation  case  is  over  1/2°K  larger 
than  for  the  corresponding  clear-sky  results.  Significant,  though  unsystematic, 
differences  in  the  rms  results  are  evident  even  into  the  stratosphere.  These 
differences  are  due  to  the  fact  that  the  lower  channels  do  contribute  signfici- 
cantly  to  the  measurement  of  temperatures  at  the  higher  levels. 

The  effects  of  moving  the  cloud  and  rain  layer  to  4 km,  in  the  case  of 
the  60  GHz  channels,  are  shown  in  Figure  3-22.  Two  significant  effects  can 
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be  noted.  The  inversion  errors  below  the  hydrometeor  layer  become  totally 
unacceptable  with  rms  errors  in  excess  of  7°K  at  the  surface,  and  3°K  at  2 km. 
The  difference  in  rms  error  between  the  two  cases  is  unsystematic  and  large. 
These  results  might  be  due  to  the  fact  that  the  data  in  the  D -Matrix  are 
unrepresentative  of  the  conditions  simulated  by  the  evaluation  samples.  More 
likely,  they  indicate  such  severe  attenuation  in  the  lower  channel,  that  all 
temperature  information  is  lost  below  4 km  and  the  D-Matrix  inversion  is  forced 
into  unrealistic  results. 

The  corresponding  inversions  for  the  118  GHz  channels  are  shown  in  Figures 
3-23  and  3-24.  These  results  show  that  there  is  a general  increase  in  rms 
error  with  the  hydrometeor  layer  elevated  to  4 km.  There  are  certain  features 
which  appear  to  be  inconsistent  with  the  results  obtained  for  the  60  GHz 
channels.  Better  temperature  inversion  results  can  be  obtained  apparently 
with  the  118  GHz  channels  for  cloudy  sky  conditions  that  is  possible  for 
clear  sky  conditions. (Figure  3-23).  However,  it  should  be  pointed  out  that 
due  to  the  effect  of  water  vapor  absorption,  very  little  information  on 
temperature  is  really  contained  in  these  channels  for  levels  of  the  atmo- 
sphere below  4 km.  Figures  3-7  and  3-8  demonstrated  that  the  lowest  peak  of 
this  set  of  weighting  functions  in  the  presence  of  water  vapor  is  at  approxi- 
mately 5 km,  due  to  the  wing  of  the  183  GHz  water  vapor  line.  Results 
obtained  below  this  level  are  not  reliable. 

3.5  Effects  of  Droplet  Size 

Two  models  of  hydrometeors  were  used  in  the  inversion  analysis  presented 
in  Section  3.4.  The  Staelin  approximation  was  used  to  model  the  effects  of 
clouds,  while  the  ERT  approximation  was  used  to  model  precipitation  effects. 
These  approximation  models  were  developed  as  an  aid  to  the  computation  of 
the  unit  volume  absorption  by  the  hydrometeor  particles.  Their  formulation 
is  not  appropriate  for  analysis  of  the  effect  of  droplet  size  independent 
of  the  effects  of  absorption  by  the  bulk  liquid  water  content  of  the  hydro- 
meteors. As  discussed  in  Section  2,  the  Staelin  approximation  for  the 
Rayleigh  regime  provides  a formulation  of  the  unit  volume  absorption  coef- 
ficient which  is  dependent  only  on  the  mass  of  water,  and  is  independent  of 
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any  drop  size  parameter.  The  liR'l  approximation  gives  the  unit  volume  absorp- 
tion coefficient  as  an  explicit  function  of  the  rainfall  rate.  Tracing  through 
the  derivations  given  in  Section  2 of  this  approximation,  it  is  possible  to 
relate  the  volume  absorption  coefficient  to  a drop  size  parameter  through  the 
Marsha  1 1 -Palmer  distribution.  However,  it  is  not  possible  to  decouple  the 
drop  size  distribution  and,  therefore,  a drop  size  parameter,  from  the  bulk 
liquid  water  content  of  the  hydrometeor.  It  is  therefore  impossible  to  per- 
form an  analysis  of  the  effects  of  changing  drop  sizes  with  the  bulk  liquid 
water  held  constant  using  this  model. 

In  order  to  satisfy  the  requirements  for  a realistic  evaluation  of 
the  effects  of  drop  size  on  brightness  temperatures,  computation  was  made 
using  the  full  Mie  approach  and  us  ng  the  analytical  distribution  of  Deirmend- 
jian  (1964).  This  distiibution  expresses  the  total  number  density  function 
N(r)  explicitly  in  terms  of  the  r iss  density  and  a drop  size  parameter,  the 
mode  radius.  The  distribution  i . given  by: 
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Computat  i ons  were  made  of  the  brightness  temperatures,  through  a tropic-' 1 
ocean  atmosphere,  using  a hvdrometeor  layer  with  uniform  mass  density  of 
1 g/mJ  extending  from  the  surface  to  4 km.  Three  sets  of  computations  were 
made,  with  mode  radius  at  20,  100  and  500  pm.  A mode  radius  of  20  pm  is 
representative  of  the  spectrum  for  cloud  droplets,  and  a mode  radius  of  500  pm 
is  representative  of  rain.  The  100  pm  mode  radius  was  selected  to  examine 
the  transition  between  cloud  and  rain. 

Results  of  the  computations  for  the  frequency  range  between  3 and  118  GHz 
are  shown  in  figure  3-25.  As  to  be  expected,  all  three  curves  converge.  For 
the  three  droplet  sizes  used  this  occurs  near  00  GHz.  The  convergence  is 
due  to  the  fact  that  as  the  drop  size  becomes  large,  relative  to  the  wavelength, 
the  Mie  efficiency  factor  tends  to  be  a limiting  value  of  2.  The  computations 
show  that  for  realistic  models  of  rain  and  cloud  drop  size  distributions,  no 
drop  size  dependence  is  expected  for  frequencies  higher  than  60  GHz. 

The  low  end  of  the  frequency  spectrum  shows  the  convergence  of  the  curves 
for  r - 20  pm  and  r,  = 100  pm.  If  the  computations  had  been  performed  for 
frequencies  lower  than  3 GHz,  the  results  would  show  a convergence  of  all 
three  curves  at  some  lower  frequency.  This  convergence  is  predicted  by  t he 
Staelin  approximation  which  states  that  as  the  drop  size  becomes  smaller, 
relative  to  wavelength,  the  unit  volume  absorption  coefficient  (and,  therefore, 
tlie  brightness  temperature)  are  independent  of  drop  size,  but  are  a function 
of  the  water  mass.  Ixtrapo 1 at ing  the  three  curves  to  lower  frequencies,  an 
approximate  convergence  frequency  of  approximately  1 GHz  can  be  obtained 
for  the  range  of  mode  radii  between  20  and  500  pm.  Below  1 GHz,  then,  no 
drop  size  effect  is  expected. 

If  the  range  of  mode  radii  between  20  and  100  pm  is  representative 
of  limiting  conditions  for  clouds,  Figure  3-25  shows  that  a cloud  drop  size 
effect  can  be  expected  in  the  range  of  frequencies  between  3 and  60  GHz,  with 
the  maximum  effect  at  approximately  10  GHz.  For  raindrops  in  the  range  of 
100  to  500  pm,  significant  drop  size  effects  are  found  in  the  frequency  range 
between  1 and  10  GHz,  with  maximum  effects  at  approximately  3 GHz. 


4.  PROBING  FOR  CLOUD  PARAMETERS 


4.1  Previous  Results 

In  previous  analyses  of  microwave  techniques  to  probe  cloud  parameters, 
the  focus  has  been  on  the  use  of  the  22.235  GHz  line  of  water  vapor.  Multi- 
channel techniques,  with  channels  on  the  wings  of  the  line  and  at  near-line 
center,  were  developed  and  analyzed  as  a means  of  separating  out  the  effects 
of  water  vapor  from  those  of  the  liquid  water  content  of  clouds.  Typical 
of  these  studies  is  that  by  Gaut,  Reifenstein  and  Chang  (1972).  Experimentally, 
the  Nimbus  E Microwave  Spectrometer  (NEMS),  with  channels  at  22.235  and  31.40 
GHz  in  addition  to  three  oxygen  chrnnels,  has  demonstrated  that  it  is  possible 
to  obtain  integrated  cloud  liquid  water  and  integrated  water  vapor  from  these 
techniques . 

However,  very  little  work  h..s  been  done  to  extend  these  techniques  to 
obtain  additional  parameters  with  regard  to  clouds.  Of  immediate  interest 
to  the  Air  Force,  for  example,  is  the  possibility  of  discrimination  between 
glaciated  and  unglaciated  clouds,  and  in  the  case  of  glaciated  clouds  the 
height  of  the  water  cloud  top  as  well  as  the  height  of  the  entire  cloud. 

It  would  also  be  of  interest  to  be  able  to  prpfile  the  vertical  distribution 
of  liquid  water  to  better  define  the  layering  structure  of  the  clouds. 

A study  was  undertaken  to  extend  the  analysis  of  remote  cloud  probing 
by  microwave  techniques  to  include  these  parameters.  Furthermore,  the 
analysis,  in  terms  of  frequency,  was  extended  to  include  the  183  GHz  water 
vapor  line.  The  simultaneous  use  of  an  infrared  channel  in  cloud  probing 
was  also  examined.  The  results  of  this  study  are  presented  in  this  section. 

4.2  Liquid  Water  Profile 

The  study  of  liquid  water  vapor  profiling  was  very  preliminary.  The 
analysis  focused  on  the  ability  of  microwave  remote  sensors  to  provide  inte- 
grated cloud  liquid  water  in  a number  of  thick  layers  in  the  atmosphere. 

As  a test,  the  atmosphere  was  divided  into  only  three  layers  as  shown  in 
the  table  that  follows. 


Preceding  page  blank 
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Layer 

Altitude  Range 

Cloud  Regime 

1 

Surface  to  3 km 

Low  c 1 oud  s 

2 

3 km  - 6 km 

Middle  clouds 

3 

6 km  - 20  km 

High  clouds 

Inversion  analysis  was  then  performed  with  various  combinations  of  channels. 
The  results  of  this  analysis  are  shown  in  Table  4-1 

The  atmospheric  regime  selected  for  the  inversion  analysis  was  that 
representative  of  midlatitude  spring.  The  sample  size  of  the  data  used  to 
generate  the  D-Matrix  was  84  and  the  evaluation  sample  size  was  36.  live 
sets  of  channels  were  analyzed,  as  shown  below. 

Channels  located  only  in  the  60  to  118  GHz 
oxygen  complexes. 

Oxygen  channels  in  combination  with  one  water 
channel 

Oxygen  channels  in  combination  with  two  water 
channels 

Channels  which  included  no  oxygen  channels 

In  examining  the  results,  it  is  interesting  to  note  that  independent  of 

the  channel  combinations,  the  best  results  in  terms  of  the  figure  of  merit 

arc  obtained  for  the  liquid  water  in  the  lowest  layer.  This  is  the  layer 

with  the  largest  variability  in  liquid  water  (a  priori  rms  variability  of 
2 

0.026  g/cm  ).  It  is  encouraging  that  with  microwave  measurements  the  rms 
error  in  the  inferred  values  for  this  layer  is  significantly  better  than 
climatology.  The  best  overall  results,  for  all  three  layers,  are  obtained 
when  the  two  sets  of  oxygen  channels  are  used  in  conjunction  with  the  two 
water  vapor  channels  at  183  and  22  Gllz,  respectively.  The  fact  that 
acceptable  results  are  obtainable  with  the  oxygen  channels  alone  (First  Case) 


First  Case 


Second  and 
Third  Cases 


Fourth  Case 


Fifth  Case 


36 
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is  probably  due  to  the  effects  of  the  183  GHz  water  vapor  line  extending  into 
the  118  GHz  region.  The  overall  results  suggest  that  it  is  possible  to 
obtain,  in  a gross  sense,  the  layering  structure  of  cloud  fields  from  micro- 
wave  remote  sensing. 

4.3  Cloud-Type  Discrimination  and  Cloud  Top  Altitude  Determination 

The  approach  to  this  analysis  was  to  assume  that  an  infrared  channel, 
similar  to  the  window  channel  on  the  DMSr,  coincidently  views  a cloud  with 
the  microwave  radiometer.  Furthermore,  it  is  assumed  that  the  clouds  arc 
optically  thick  and  of  sufficient  horizontal  extent  such  that  the  effective 
blackbody  temperature,  T^,  measured  by  the  infrared  channel,  approximates 
for  all  practical  purposes  the  actual  physical  temperature  of  the  air  at 
equilibrium  with  the  cloud  top.  This  cloud  top  temperature  would  correspond 
to  the  temperature  of  the  physical  top  of  the  cloud,  regardless  of  the  phase 
characteristics  of  the  cloud  top  layer.  If  a microwave  frequency  is  found 
such  that  it  is  totally  transparent  to  ice,  but  highly  opaque  to  liquid 
water,  its  brightness  temperature  would  be  well  correlated  with  the  top  of 
the  liquid  portion  of  the  cloud.  The  difference  between  the  effective  black- 
body  temperature  measured  by  the  IR  channel  and  the  brightness  temperature 
measured  by  this  microwave  channel  would  provide  information  on  the  presence 
of  a glaciated  layer.  Furthermore,  the  measured  temperatures  themselves 
would  provide  data  on  the  heights  of  the  water  cloud  top  and  the  ice  cloud 
top. 

Analysis  was  performed  to  find  a suitable  microwave  channel  to  work 
with  the  infrared  channel.  The  analysis  focused  on  the  use  of  the  183  GHz 
water  vapor  line  which  is  much  more  opaque  than  the  22  GHz  line.  Furthermore 
computations  shown  in  Figure  3-25  demonstrate  that  at  the  higher  frequencies 
the  brightness  temperature  for  a hydrometeor  layer  is  independent  of  drop 
size.  Computations  were  made  for  three  frequencies  at  140  and  170  GHz  on 
the  wings  of  the  line  and  at  183  GHz  near  line  center.  For  the  sake  of 
analysis,  three  cloud  model  conditions  were  assumed:  a 4 km  cloud  layer 

extending  to  the  surface,  the  same  cloud  layer  but  with  a 2 km  layer  of  ice 
-ibove  it,  and  an  8 km  water  cloud  with  a 4 km  layer  of  ice.  The  temperatures 
of  the  cloud  tops  for  these  models  are  summarized  in  Table  4-2. 
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TABU:  4-2 

SUMMARY  OF  THE  TEMPERATURE  CHARACTERISTICS  OF  CLOUD  MODELS 

USED  FOR  ANALYSIS 


Cloud  Me  'el 
Cloud  Model 

Temperature  v k) 

Liqu;d 

I cc 

4 km  water 

275.6 

NA 

275.  (. 

4 km  and  2 km  ice 

275.6 

263.7 

263.7 

8 km  water  and  4 km  ice 

250.0 

223.5 

223.3 

Results  of  these  computations  are  shown  in  Figures  4-la  and  4-lb. 
Figure  4-la  shows  that  the  brightness  temperature  at  the  selected  fre- 
quencies are  only  weakly  dependent  on  cloud  liquid  water  content.  Fur- 
thermore, the  brightness  temperature  at  170  GHz  is  almost  identical  to 
the  physical  temperature  of  the  cloud  top.  Similar  computations  for  an 
8 km  liquid  water  cloud,  with  cloud  top  at  750°K,  also  showed  that  the 
brightness  temperature  at  170  GHz  corresponds  quite  well  to  the  physical 
temperature  of  the  cloud  top.  Further  analysis  is  required  to  demon- 
strate whether  or  not  this  correspondence  holds  in  the  general  ease. 

The  effects  of  an  ice  layer  on  IR  effective  temperature  and  micro- 
wave brightness  temperature  is  shown  in  Figure  4-lb.  For  these  calcu- 
lations, the  liquid  cloud  properties  were  held  constant  with  a liquid 

3 

cloud  density  of  0.6  gm/cm  (corresponding  to  an  integrated  liquid  water 

2 

content  of  0.25  gm/cm“  for  the  4 km  cloud  layer).  The  ice  density  in 
the  2 km  ice  layer  was  changed;  thus  effectively  changing  the  equivalent 
liquid  water  content  of  the  ice  layer  without  changing  the  geometric 
properties  of  the  cloud.  The  plotted  1R  effective  cloud  top  tempera- 
ture again  assumed  that  the  ice  layer  is  essentially  opaque  in  the 
IR;  thus  the  effective  blackbody  temperature  measured  in  the  1R  corres- 
ponds to  the  geometric  top  of  the  entire  cloud. 

In  comparing  the  results  shown  in  the  two  figures,  it  is  seen  that 
within  the  limits  of  the  models  used  in  the  calculations,  the  microwave 
brightness  temperature  of  clouds  at  these  selected  frequencies  are  not 
strongly  dependent  on  the  density  of  an  ice  layer.  These  results  imply 
that  brightness  temperatures  obtained  using  channels  at  these  high 
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frequencies,  when  used  with  concurrent  IR  blackbody  temperatures,  can  be 
of  use  to: 

a)  differentiate  between  glaciated  from  unglaciated  clouds,  and 

b)  estimate  the  height  of  both  the  top  of  the  cloud  and  the  top 
of  the  liquid  cloud  in  the  presence  of  an  ice  layer. 

Computations  using  an  8 km  water  cloud  with  and  without  a 4 km  ice  layer 
support  the  results  shown  here  for  the  4 km  cloud  case.  Further  research, 
however,  is  needed  to  pinpoint  a useful  operational  technique  applicable 
to  clouds  with  much  more  complex  structure  than  those  modelled  in  this 
study . 
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Brightness  Temperature  (°K) 


102*0 


1'itegrated  Liquid  Water  Content  (g  cm-2) 

a)  Liquid  Clouds  with  Cloud  Top  Temperature 
at  275.6°  K 


0 0.1  0.2  0.3  0.4  0.5  0.6 


Ice  Density  (g  cm-3) 

b)  Glaciated  Clouds  with  Liquid  Cloud  Top 
Temperature  at  275. 6°K  and  Ice  Cloud 
Top  Temperature  at  263.7  K 


Figure  4-1  Microwave  and  Infrared  Measurements  of  Glaciated  and  Unglaciated  Clouds 
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5.  SUMMARY  AND  CONCLUSIONS 


An  approximate  model  for  the  computation  of  the  unit  extinction  coeffici- 
ent for  rain  has  been  developed  to  facilitate  simulation  and  inversion  analyses. 
In  terms  of  the  interaction,  the  model  assumes  that  the  efficiency  factor  as 
a function  of  drop  radius  can  be  approximated  by  three  regimes:  Rayleigh 

expression,  transition,  and  large  sphere.  Analytical  expressions  have  been 
developed  for  these  regions  such  that  the  convolution  of  the  efficiency  factor 
with  the  Marshall-Palmer  distribution,  by  which  the  drop  size  distribution  is 
expressed,  can  be  integrated  in  closed  form.  The  use  of  the  Marshall -Palmer 
distribution,  taken  from  meteorological  radar  studies,  results  in  an  expres- 
sion which  relates  the  synoptic  parameter  rainfall  rate  to  the  precipitation 
extinction  in  a given  layer.  In  addition,  temperature  dependence  and  wave- 
length dependence  are  analytically  expressed  providing  a realistic  layer 
parameterization  which  is  amenable  to  numerical  computation.  It  has  been 
shown  that  the  computed  extinction  coefficients  faithfully  reproduce  full 
Mie  computations  within  a reasonable  accuracy.  This  development  enables 
computations  of  the  extinction  effect  of  precipitation  in  atmospheres  with 
realistic  temperature  structure  without  resorting  to  the  full  Mie  expansion. 

This  model  for  rain,  and  the  Staelin  model  for  clouds,  were  used  to 
investigate  the  effects  of  hydrometeors  on  temperature  inversion  accuracies. 

The  inversion  algorithm  used  was  the  statistical  method,  and  it  was  applied 
to  two  sets  of  channels:  one  in  the  60  GHz  complex  and  the  other  in  the  118 

GHz  complex  of  oxygen.  In  general,  the  results  show  that  clouds  and  rain  do 
have  significant  deteriorating  effects  on  temperature  inversion  accuracies. 
These  effects  are  much  more  pronounced  in  the  case  of  the  channels  in  the 
118  GHz  region.  Furthermore,  the  hydrometeors  tend  to  produce  errors  in 
inversion  at  altitudes  many  kilometers  above  the  tops  of  the  layers. 

The  effects  of  drop  size  on  the  brightness  temperature  of  hydrometeor 
layers  were  separately  analyzed  using  the  full  Mie  expansion  and  the 
Diermendjian  expression  for  drop  size  distribution.  For  the  model  cloud 
analyzed,  the  results  indicate  that  the  effect  of  drop  size  on  brightness 
temperature  is  separable  into  two  regimes:  rain  and  clouds.  For  rain,  the 

effects  of  drop  size  are  confined  to  the  frequency  range  between  1 and  10 
GHz,  with  maximum  effects  at  3 GHz.  For  cloud  droplets,  the  effects  of 
drop  size  can  be  noted  in  the  range  of  frequencies  from  3 to  60  GHz,  with 
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maximum  effect  at  10  GHz.  It  would  appear  that  the  frequency  region  best 
suited  for  remote  probing  of  drop  size  parameters  is  found  between  3 and  10 
GHz,  provided  that  some  independent  measurements  can  be  made  to  separate  out 
the  effects  of  bulk  liquid  water  of  the  hydrometeors. 

Inversion  analysis  performed  in  this  study  suggests  that  with  the  right 
combination  of  channels,  it  is  possible  to  infer  the  layered  structure  of 
clouds.  Furthermore,  by  using  a very  opaque  channel  located  in  the  wing  of 
the  183  GHz  water  vapor  line  in  conjunction  with  an  infrared  channel,  the 
possibility  to  differentiate  glaciated  clouds  from  water  clouds  exists. 

The  same  channels  can  provide  additional  information  about  the  height  of 
the  top  of  the  liquid  cloud  as  well  as  that  of  the  ice  layer. 
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